THIRD-ORDER NONLINEAR DISPERSION PDES: 
SHOCKS, RAREFACTION, AND BLOW-UP WAVES 



VICTOR A. GALAKTIONOV AND STANISLAV I. POHOZAEV 

Abstract. Formation of shocks and collapse of other discontinuities via smooth rar- 
efaction waves, as well as blow-up phenomena, for the third-order nonlinear dispersion 
PDE, as the key model, 

(0.1) ut = {uu^)^^ in R X R+, 

are studied. Two basic Riemann's problems for (jO.ip with initial data 

^Tix) = Tsignx, 

are shown to create the shock {u{x,t) = S-{x)) and smooth rarefaction (for data S+) 
waves, respectively. To this end, various blow-up and global similarity solutions of (jO.ip 
revealing the fine structure of shock and rarefaction profiles are constructed. Next, 
eigenfunction and nonlinear capacity techniques for proving blow-up are developed. 
The study of (jO.ip reveals similarities with entropy theory for scalar conservation laws, 

Ut + UUx — 0, 

developed by Oleinik and Kruzhkov (equations for x G M^) in the 1950-60s. 



Dedicated to the memory of Professors O.A. Oleinik and S.N. Kruzhkov 

This is the first paper on shock wave theory for higher-order nonlinear dispersion equa- 
tions, which will be continued in three more ones [121 [IS] and 



1. Introduction: nonlinear dispersion PDEs and main directions of study 

1.1. NDEs: nonlinear dispersion equations in application and general PDE 
theory. In the present paper, we develop basic aspects of singularity and existence- 
uniqueness theory for odd-order nonlinear dispersion (or dispersive) PDEs, which, for 
short, we call the NDEs. The canonical for us model is the third- order quadratic NDE 
(the NDE-3) 

(1.1) Ut = A{u) = {uux)xx = uuxxx + ^UxUxx in M x (0, T), T > 0. 
We pose for (11.11) the Cauchy problem with locally integrable initial data 

(1.2) m(x,0) = uo(a;) in M. 
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Often we assume that Uq is bounded and/or compactly supported. We will also deal with 
the initial-boundary values problem in {—L, L) x with Dirichlet boundary conditions. 

On integrable NDEs related to water wave theory. Concerning applications of equations 
such as (11.11) . it is customary that various odd-order PDEs appear in classic theory of 
integrable PDEs, such as the KdV equation., 

(1.3) ut + uu^ = u^^^, 
and the fifth-order KdV equation, 

Ut + Uxxxxx + 30 U^Ux + 20 UxUxx + 10 UUxxx = 0, 

and others from shallow water theory. The quasilinear Harry Dym equation 

(1.4) Ut = U^Uxxx, 

which also belongs to the NDE family, is one of the most exotic integrable soliton equa- 
tions; see [TTl § 4.7] for survey and references therein. Integrable equation theory produced 
various hierarchies of quasilinear higher-order NDEs, such as the fifth-order Kawamoto 
equation [22j 

Uf U Uxxxxx ~\~ 5 U UxUxxxx ~l~ 10 It UxxUxxx- 

Quasilinear integrable extensions are admitted by Lax's seventh-order KdV equation 

Ut + {35U^ + 70[u^Uxx + u{Uxy] + 7[2uUxxxx + ^iUxxY + ^UxUxxx] + Uxxxxxx}^ = 0, 

and by the seventh-order Sawada-Kotara equation 

Ut + {63m^ + 63[2u\xx + u{uxf] + 21 [uUxxxx + (uxx)'^ + UxUxxx] + Uxxxxxx}^ = 0; 

see references in [ITJ p. 234]. We will continue this short survey on existence-uniqueness 
for some integrable odd-order PDEs in Section [T31 

Compact pattern formation and NDEs. Returning to lowest third-order NDEs that are 
not integrable, we will also study the Rosenau-Hyman (RH) equation 

(1-5) Ut = {U^)xxx + {u'^)x 

which has special important applications as a widely used model of the effects of nonlinear 
dispersion in the pattern formation in liquid drops [34] . It is the 7^^(2,2) equation from 
the general K{m, n) family of the following NDEs: 

(1.6) Ut = {u'')xxx + {u"')x {u>0), 

that also describe various phenomena of compact pattern formation for m, n > 1, [301131]. 
Such PDEs also appear in curve motion and shortening flows [53]. Similar to well-known 
parabolic models of porous medium type, the K{m, n) equation (11.61) with > 1 is 
degenerated at u = 0, and therefore may exhibit finite speed of propagation and admit 
solutions with finite interfaces. The crucial advantage of the RH equation (II. 5p is that it 
possesses explicit moving compactly supported soliton-type solutions, called compactons 
[M] , which are travelling wave (TW) solutions. We will check entropy properties of 
compactons for various NDEs. 
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Various families of quasilinear third-order KdV-type equations can be found in [3], 
where further references concerning such PDEs and their exact solutions are given. Higher- 
order generalized KdV equations are of increasing interest; see e.g., the quintic KdV equa- 
tion in [20] and ^41j, where the seventh-order PDEs are studied. For the K{2, 2) equation 
(11. 5p . the compacton solutions were constructed in [30]. More general B{m, k) equations 
(indeed, coinciding with the K{m, k) after scaling) 

also admit simple semi-compacton solutions [32], as well as the Kq{m,u) nonlinear dis- 
persion equation (another nonlinear extension of the KdV) pD] 

Ut + {u^')x + [u'-'^iu'^u,),], = 0. 

Setting m = 2 and uj = ^ yields a typical quadratic PDE 

Ut + {u'^)x + uu^^^ + 2m^m^^ = 

possessing solutions on standard trigonometric-exponential subspaces, where 

u{x, t) = Co(t) + Ci{t) cos Xx + C2{t) sin \x 

and {Co,Ci,C2} solve a 3D nonlinear dynamical system. Combining the K{m,n) and 
B{m,k) equations gives the dispersive-dissipativity entity DD{k,m,n) [32] 

Ut + ci(m™)x + {u'')xxx = Ku'')xx 

that can also admit solutions on invariant subspaces for some values of parameters. 
For the fifth-order NDEs, such as 

(1.7) Ut = a{u^)xxccxx + P{u'^)xxx + l{u'^)x in M X IR+, 

compacton solutions were first constructed in [Sj, where the more general K{m,n,p) 
family of PDEs 

Ut + Piiu^x + f32{unxxx + PsDliuP) =0 {m,n,p> 1) 

was introduced. Some of these equations will be treated later on. Equation (11.71) is 
also associated with the family Q{l,m,n) of more general quintic evolution PDEs with 
nonlinear dispersion, 

(1.8) Ut + a(u-+i), + u [u{u^),,] ^ + 6 [m(u%...] ^ = 0, 

possessing multi-hump, compact solitary solutions [36j- 

Concerning higher-order in time quasilinear PDEs, let us mention a generalization of 
the combined dissipative double- dispersive (CDDD) equation (see, e.g., [28]) 

(1-9) utt = au^^xx + (3u^xtt + l{u^)xxxxt + S{u^)xxt + £{u^)t, 

and also the nonlinear modified dispersive Klein-Gordon equation {mKG{l,n, k)), 

(1.10) + a(M")^.^. + 6(m^)^^^^ = 0, n,k>l (m > 0); 
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see some exact TW solutions in [2T]. For 6 > 0, (11.101) is of hyperbolic (or Boussinesq) 
type in the class of nonnegative solutions. Let us also mention a related family of 2D 
dispersive Boussinesq equations denoted by B{m,n, k,p) [ID], 



{u"')tt + + (3iu^)xxxx + liu^)yyyy = in X R. 

See p!7t Ch. 4-6] for more references and examples of exact solutions on invariant subspaces 
of NDEs of various types and orders. 

NDEs in general PDE theory. In the framework of general theory of nonlinear evolution 
PDEs of the first order in time, the NDE (II. ip appears the third in the following ordered 
list of canonical evolution quasilinear degenerate equations: 

(1.11) Ut = — I {u^)x (the conservation law), 

(1.12) Ut = \ {v?)xx (the porous medium equation), 

(1.13) Ut = \ {u'^)xxx (the nonlinear dispersion equation), 

(1.14) Ut = — \ (\u\u)xxxx (the 4*'^-order nonlinear diffusion equation). 

In (11.141) . the quadratic nonlinearity v? is replaced by the monotone one \u\u in order to 
keep the parabolicity on solutions of changing sign. The same can be done in the PME 
(11.121) . though this classic version is parabolic on nonnegative solutions, a property that 
is preserved by the Maximum Principle. The further extension of the list by including 

(1.15) ut = -\{u')xxxxx (theNDE-5) and 

(1.16) '^t = \ {\u\u)xxxxxx (the 6*'^-order nonlinear diffusion equation) 

is not that essential since these PDEs belong to the same families as (I1.13P and (11.140 
respectively with similar covering mathematical concepts (but indeed more difficult). 

Mathematical theory of first two equations, (II. lip (see detailed survey and references 
below) and (I1.12p (for quoting main achievements of PME theory developed in the 1950- 
80s, see e.g., [18], Ch. 2]), was essentially completed in the twentieth century. It is curious 
that looking more difficult the fourth-order nonlinear diffusion equation (I1.14p has a mono- 
tone operator in so the Cauchy problem admits a unique weak solution as follows 
from classic theory; see Lions [231 Ch. 2]. Of course, some other qualitative properties of 
solutions of (I1.14P are more difficult and remain open still. 

It turns out that, rather surprisingly, the third order NDE (ll.lSp is the only one in this 
basic list that has rather obscure understanding and lacking of a reliable mathematical 
basis concerning generic singularities, shocks, rarefaction waves, and entropy-like theory. 
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1.2. Mathematical preliminaries: analogies with conservation laws, Riemann's 
problems for basic shocks S^{x) and if (=fx), and first results. As a key feature of 
our analysis, equation (11. ip inherits clear similarities of the behaviour for the first-order 
conservation laws such as Euler's equation (same as (11.111) ) from gas dynamics 

(1.17) ut + uu^ = in M X M+, 

whose entropy theory was created by Oleinik j26l [27] and Kruzhkov ^23j (equations in 
M^) in the 1950-60s; see details on the history, main results, and modern developments 
in the well-known monog raphs [I El EHfl. 

As for (I1.17p . in view of the full divergence of the equation (II. ip . it is natural to 
define weak solutions. For convenience, we present here a standard definition mentioning 
that, in fact, the concept of weak solutions for NDEs even in fully divergent form is not 
entirely consistent, to say nothing about other non-divergent equations admitting no weak 
formulation at all. 

Definition [T].l. A function u = u{x,t) is a weak solution of (II. ip . (I1.2p if 

(i) «2eLL(Mx(o,r)), 

(ii) u satisfies (II. ip in the weak sense: for any test function ip{x,t) G C^(]R x (0,T)), 

(1-18) JJ Ulft = I JJ U'^fxxx, 

and (iii) satisfies the initial condition (ll.2p in the sense of distributions, 

(1.19) ess limf_,o / u{x,t)ip{x) = J Uo{x)iIj{x) for any ip G C^(M). 

We will show that, for some data, (II. 2p is also true with convergence in Lj'Q(,(R). This 
topology is natural for conservation laws (I1.17p . but for the NDE (II. ip such a convergence 
is not straightforward and demands a precise knowing of the structure of oscillatory "tails" 
of similarity rarefaction solutions to be studied in detail. The assumption T < oo is often 
essential, since, unlike (I1.17p . the NDE (II. ID can produce complete blow-up from bounded 
data. 

Thus, again similar to (I1.17p . one observes a typical difficulty: according to Definition 
[Hi, both discontinuous step-like functions 

{1 for a; > 0, 
for X = 0, 
— 1 for X < 0, 

are weak stationary solutions of (ll.ip satisfying 
(1.21) iu^)xxx = 

in the weak sense, since in (11.180 m^(x) = S^{x) = 1 is C'^ smooth (and analytic). 
Again referring to entropy theory for conservation laws (I1.17p . it is well-known that 

^ S_{x) is the entropy shock wave, and 

S+{x) is not an entropy solution. 



First results on the formation of shocks were performed by Riemann in 1858 [52]; see [3]. 
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This means that 



(1.23) M„(x, t) = S^{x) = — signx 

is the unique entropy solution of the PDE fll.l7p with the same initial data S-{x). On 
the contrary, taking 5+ initial data yields the rarefaction wave with a simple similarity 
piece-wise linear structure 

{—1 for X < —t, 
J for |a;| < t, 
1 for X > t. 

First aims of singularity theory. Thus, as in classic conservation law theory, the above 
discussion formulate our first two aims: 

• Our first one is to justify the same classification of main two Riemann's problems 
with data fll.20p for the NDE (11. II) and to construct the corresponding rarefaction wave 
for Sj^{x), as an analogy of (11.241) for the conservation law; 

• The second aim is to describe evolution blow-up formation for (II. ip of another shock 
that is the reflected Heaviside function given by 

1 for X < 0, 
for X > 



;i.25) H{ 



(iJ(0) = 0.4197... will be determined by construction). It turns out that this third Rie- 
mann problem for (11.11) is also solved in a blow-up similarity fashion with the solution 
exhibiting a. finite interface. As a crucial feature, note that H{—x) is not a weak stationary 
solution of (II. ip . In general, we therefore would like to claim the following rule: 

(1.26) shock waves of higher-order NDEs are not necessarily weak solutions. 

Of course, this is obvious for non-divergent NDEs or those that are fully nonlinear (we 
will consider such equations), but this also happens for fully divergent PDEs such as (II. ip . 

1.3. NDEs are not hyperbolic systems. It should be noted that by no means (II. ip 
can be treated as a hyperbolic system, for which fully developed theory in ID is now 
available; see e.g., Bressan [2]. Indeed, rewriting (II. ip as a first-order system yields 

(1.27) 

Evidently, it is not a hyperbolic system that represents a combination of a single first- 
order evolution PDE with two stationary equations. In a natural formal sense, (I1.27P can 
be considered as the limit as e — of a more standard first-order system such as 

= w^, 

1.28) ivt=^-(v-uu 




e 

Wt = \ {w — V, 
6 



Nevertheless, the principal differential operator in fll.28p is not hyperbolic for any |e| <^ 1 
that is easily seen by calculating eigenvalues of the matrix, 

-A 1 

'1.29) - A 

-i -A 



^ A^ = 4. 



The algebraic equation in (11.291) admits complex roots for any m 7^ that reflects the 
actual highly oscillatory properties of solutions of both systems (11.291) and (11.281) and 
the original PDE (II. ip . Note that the generic behaviour of small solutions of the KdV 
equation (II. 3p with linear dispersion mechanism is also oscillatory. This is associated with 
the fundamental solution of the corresponding operator Dt — given by Airy function 
of changing sign, which is not absolutely integrable on the whole line x G M for any t > 0. 

1.4. Conservation law in H^^. On the other hand, equation fll.ip can be written as 

(1.30) {-Dl)-\t + uu^ = 0, 

with the following usual definition of the linear operator (— Z^^.)^^ > in R: 



-00 



{-Dlr^v = g, if g" = -v, gi±c^) = {Jv = 0), i.e., 
. . —00 
1.31 

00 y X 

g{x) = J J v{z)dzdy, provided that v e L'^{R), J v e L^{R). 

X —00 —00 

In a bounded interval [—L, L] to be also used, the definition is standard, with g{^L) = 0. 

The form (ll.30p makes it possible to get the first a priori uniform bound on solutions 
for data Uo E H~^{R): multiplying fOD]l by u in we get the conservation law 
(1.32) 1 d ||^^^)||2^_^ ^ ||^(i)||2^_^^^^^||^^||2^_^ forallt>0. 

Notice that this is weaker than the estimate for the conservation law (I1.17p . for which 
(I1.32P takes place in the topology of L^(R). 

It seems that the representation fll.30p and the estimate (11.320 do not help neither better 
understanding of formation of singularities nor entropy essence of solutions. Nevertheless, 
(ll.30p convinces to classify the NDE (11.10 as a non-local conservation law posed in H~^. 

1.5. Plan of the paper: known results and main directions of our singular- 
ity study. We mention again that, surprisingly, quite a little is known about suitable 
mathematics of shocks, rarefaction waves, and entropy essence of the NDEs such as (II. ip . 
Moreover, even for sufficiently smooth, continuous solutions, it is not exaggeration to say 
that, for (II. ip and the compacton equation (11.50 . even basic facts concerning proper pos- 
ing the Cauchy problem remained obscure, especially, existence of discontinuous shodO 



^ It is relevant commenting that a few years ago, one of the first author's paper on shock waves for 
the NDE submitted to a respectable journal partially dealing with nonlinear PDEs received the 

accusation that "... this equation cannot admit discontinuous solutions at all as the heat equation". It 
seems the Referee was confused observing in (11.11) too many (three) spatial derivatives, while the heat 
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and rarefaction waves, as well as compactons as compactly supported solutions (in the 
twenty-first century!). 

In the present paper, we describe various types of shock, rarefaction, and blow-up 
waves for NDEs under consideration. In the next paper |12j, we are going to propose 
key concepts for developing adequate mathematics of NDEs with shocks, which will be 
concluded by revealing connections with other classes of nonlinear degenerate PDEs. It 
turns out that some NDE concepts has definite reliable common roots and can be put 
into the framework of much better developed theory of quasilinear parabolic equations. 

Back to integrable models: existence and uniqueness. In fact, modern mathematical the- 
ory of odd-order quasilinear PDEs is originated and continues to be strongly connected 
with the class of integrable equations. Special advantages of integrability by using the 
inverse scattering transform method. Lax pairs, Liouville transformations, and other ex- 
plicit algebraic manipulations made it possible to create rather complete theory for some 
of these difficult quasilinear PDEs. Nowadays, well-developed theory and most of rig- 
orous results on existence, uniqueness, and various singularity and non-differentiability 
properties are associated with NDE-type integrable models such as Fuchssteiner-Fokas- 
Camassa-Holm (FFCH) equation 

(1.33) (/ - Dl)ut = -3uux + lu^Uxx + tt^xxx = -{I - Dl){uu^) - [u^ + \ [uxY] ■ 

Equation (11.331) is an asymptotic model describing the wave dynamics at the free surface 
of fluids under gravity. It is derived from Euler equations for inviscid fluids under the 
long wave asymptotics of shallow water behavior (where the function u is the height of 
the water above a fiat bottom). Applying to ( 11.33P the integral operator (J — D^)~^ with 
the L^-kernel uj{s) = ^e"''^' > 0, reduces it, for a class of solutions, to the conservation 
law (II. ip with a compact /irsi- order perturbation, 

(1.34) ut + uu, = -[uj*{u' + \{u^f)]^. 

Almost all mathematical results (including entropy inequalities and Oleinik's condition 
(E)) have been obtained by using this integral representation of the FFCM equation; see 
a long list of references given in [TTl p. 232]. 

There is another integrable PDE from the family with third-order quadratic operators, 

(1.35) Ut - u^xt = auux + PuxUxx + uuxxx (a, /5 e M), 

where a = —3 and (3 = 2 yields the FFCH equation (ll.33p . This is the Degasperis-Procesi 
equation for a = — 4 and P = 3, 

(1.36) Ut - u^^t = -^uu^ + 3m^m^.^. + uu^^^. 

On existence, uniqueness (of entropy solutions in fl BV), parabolic e-regularization, 
Oleinik's entropy estimate, and generalized PDEs, see [5J. Besides (11.331) and (11.361) . the 
family (I1.35P does not contain other integrable entries. A list of more applied papers 
related to various NDEs is also available in [17[ Ch. 4]. 

equation ut = Uxx has two only. Of course, any odd, 2m + 1 (say, 2007 for m — 1003), number of 
derivatives in ut = {—l)'"^^^D'jj"'{uUx) will lead to discontinuous shocks for some initial data. 
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Main directions of study. Concerning the simple canonical model (11.11) . which is not inte- 
grable and does not admit a reduction like fll.34p with a compact first-order perturbation, 
we do the following research: 

• (i) To check whether (11. ip admits discontinuous shock wave solutions, and which 
smooth similarity solutions can create shocks in finite time. 

• (ii) To verify which discontinuous solutions are entropy in a natural sense, and which 
are not that give rise to smooth rarefaction waves. To this end, we apply the idea of 
smooth deformations of shocks to see which ones are stable. In particular, we show that 
two basic Riemann's problems for (11.11) with initial data (11.201) correspond the shock 

and rarefaction (S'+) waves respectively. This coincide with the classification for the 
conservation law (I1.17p . We also solve the third Riemann problem with the Heaviside 
function (11.250 that turns out to be entropy. Similarly, the reflected data H{x) are shown 
to produce the corresponding smooth rarefaction wave. 

• (iii) For these purposes, construct various self-similar solutions of (II. ip describing for- 
mation of shocks and collapse of non-evolution discontinuities (i.e., leading to rarefaction 
waves); here we also use Gel'fand's concept (1959) [19, §§ 2,8] of G-admissible solutions 
of higher-order ODEs. 

• (iv) As a consequence, to prove that (11. ip describes processes with finite propagation 
of perturbations (indeed, this was well-known for decades but, it seems, was not suitably 
treated mathematically) . 

In the forthcoming paper [12], we continue our study and: 

• (v) Develop local existence and uniqueness theory for the NDE (II. ip and introduce 
the concept of the 5-entropy test. 

• (vi) For the RH equation such as (11.50 . we prove that Rosenau's compacton solutions 
are both entropy and G-admissible. 

Some of related questions and results were previously discussed in a more applied and 
formal fashion in [Til § 7] and [171 Ch. 4], where another idea of regular approximations 
was under scrutiny. Namely, proper solutions therein were obtained in the limit e ^ 0^ 
of classical solutions {^e} of a family of uniformly parabolic equations with same data 

(1.37) ut = {uux)xx - euxxxx in M x M+ (e: > 0). 

In particular, it was shown that a direct verification that the e-approximation (11.370 yields 
as e — > the correct Kruzhkov's-type entropy solution leads to difficult open problems. 
Further extensions to other three- and higher-order NDEs can be found in [131 [12] . 

2. First blow-up results by two methods 

In this section, we begin our study of blow-up singularities that can be generated by 
the NDE (ll.ip and others. We pose an IBVP and demonstrate two different methods 
to prove global nonexistence of a class of solutions. Note that these results differ from 
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those that have been known in the theory of more classic semilinear dispersive equations 
(SDEs) such as 

(2.1) iut + R{-iD,)u-Q{-iD.,){uP+') = 0, p>l, 

where, typically, the functions P and P are polynomials in one variables with real coeffi- 
cients; see m for the results and a survey. 

2.1. A modified eigenfunction method. We begin by deriving simple ordinary dif- 
ferential inequalities (ODIs) for Fourier-like expansion coefficients of solutions. Here 
we perform a convenient adaptation of Kaplan's eigenfunction method (1963) to odd, 
(2m-|- l)th-order nonlinear PDEs. This completes our analysis performed for even, 2mth- 
order parabolic equations in [15^ § 5], where key references and other basic results can 
be found. Another alternative approach to nonexistence for such PDEs is developed in 
author's resent paper [16]. Later on, in Section HI we will show that looking analogous 
L^-blow-up is actually possible and can be driven by self-similar solutions. 

For convenience and simplicity, we consider (11. ip for x < and take initial data 

(2.2) Uo{x)<0 (no(x)^O). 
We also assume the following conditions at the origin: 

(2.3) M = 0, (u^)^ > and the flux (u^)^^ < at x = 0. 

Note that all these conditions hold for the Cauchy problem if the right-hand interface of 
the given compactly supported solution belong to {x < 0} for t G (0, T). We will describe 
such solutions with finite interfaces later on. On the other hand, we may pose an initial- 
boundary value problem (IBVP) for (11. ip in M_ x (0, T) with the necessary conditions 
at X = such that (12.31) holds. Not discussing here the local well-posedness of such an 
IBVP, we will show that some of its solutions will then blow-up in finite time. 
We next fix a constant L > and choose the "eigenfunction" as follows: 

(2.4) 0(x) = -{x + Lf < for xE (0, L), 0(-L) = 0'(-L) = 0"(-L) = 0. 
Here satisfies the obvious "eigenfunction" equation 

(2.5) P30 = (p"'{x) = -6 on (-L, 0) (0 < 0), 

plus the above boundary conditions. Here the linear operator P3 = is naturally 
associated with the NDE (II. ip . Actually, any satisfying the differential inequality 

(2.6) <P"'{x)<-Co (0<O), 

where Co > is a constant, would successfully fit into our further analysis. 
Let us introduce the corresponding expansion coefficient of the solution 


J{t) = J u{x,t)(p{x) dx = ~ J u{x,t){x + LY dx, 

(2.7) „ 

where J(0) = — J uo{x){x + L)^ dx > 0. 
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Multiplying (11.11) by <f){x) and integrating over {—L, 0) by parts three times, by using the 
properties in (12.31) . (12.41) . and (12.51) . one obtains that 



J'{t) = i - + m2 0^^] _^ _ 1 J ^2(^^^ t)(j)"'{x) dx 

(2-8) 

> -| / M2(a;^t)0'"(a;)da; = 3 / u\x,t)dx. 

-L -L 

We next use the Holder inequality to estimate the right-hand side in (12.81) : 





3 Ju\x,t)dx>3 J u\x,t)^-^dx 

-L -L 

f {\x + Lf\u{x,t)\)^ldx > ^ ( / \x + L\^\u{x,t)\dxY > jjr-P{t). 



-L -L 

Thus, we obtain the following simple quadratic ODI for the expansion coefficient 

(2.9) J\t) > jjfJ'^it) for t > 0, J(0) = Jo > 0. 
Integration over (0, t) yields a typical blow-up, 

(2.10) J{t) > 3(^^ with the blow-up time T <To = j^. 

Recall again that this blow-up can happen if the flux {u'^)xx at x = satisfles the sign 
condition in (12. 3p . Otherwise, if it is violated and the flux gets positive, this can prevent 
any blow-up in (but cannot prevent a "weaker" gradient blow-up of Ux to be studied 
in greater detail). 

This approach to blow-up admits natural adaptation to flfth- and higher-order NDEs 

(2.11) Uf (yUUxx^x = '^'^xxx ~l~ '^x'^xxi 

which we call the NDE-(2,1), where 2 and 1 stand for the number of the internal and 
external derivatives in this differential form. For second- and third-order in time NDEs 

(2.12) utt = {uux)xx and Um = {uUx)xx- 
the same integration by parts leads, instead of (12.91) . to ODIs 

(2.13) J"{t) > ^J\t) and J"'{t) > ^J\t). 

Here blow-up depends on initial values Jo = J(0), J'(0), and J"(0). Finally, these formally 
deflne respectively the following blow-up rates: 

(2.14) j(t)^ 1 and J{t) 
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2.2. Nonlinear capacity method: another version of blow-up. We consider the 
second-order in time NDE in fl2.12l) in IR+ x , 

(2.15) Utt = \ {V?)xxx 

and, for simplicity, assume all three zero conditions (12. 3p at the origin. We now choose a 
different cut-off function, 

0(a;,t) = 0o(t)(^ - x)^ where 0o > 0, 0o(O) = 0[,(O) = 0, 0o(T) = 0[,(T) = 0. 

Integrating by parts yields 

L L 

I / {u^)xxx4>{x, t)dx = 3 f u^(po(t) dx, 


T L L T L 

J J Uu(j)dxdt = — J ut{x,0){x — L)^ dx + J J u(j)Q{t){x — L)^ dx dt. 



Hence, we obtain the integral identity 

T L L T L 

(2.16) 3/ ^ u^(t)Q{t)dxdt = - ^ Ut{x,Q){x- Lfdx + J J u(t)'^{t){x - Lf dxdt. 





By Young's inequality, setting for convenience 0i(x) = (x — L)^, 

t<|//-Vo + i//^,„ 



so that fl2.16p implies 

(2.17) I / / ^Vo < J/ - Jutix, 0)Mx) dx. 

DO 

We next perform scaling 

t I— > r = |t and replace (poit) h-> 0o(t), and then 
/ / (^) d-dt = / ^ dt / 0, (a:) dx 



- - J ^^t^r=± - — T^^O, Co - J dr < OO. 

Finally, we obtain the estimate 

T L L 

(2.18) |//?^Vodxdt < 2^ - Jo, where Jq = ^ Ut{x,Q){L - xf dx. 



Thus, if Jo > 0, the solution must blow-up in finite time Tq satisfying To < (^;|^) ^. It is 
curious that this blow-up conclusion depends on second data Ut{x, 0) and does not involve 
the first initial function m(x, 0). 
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3. Shock and rarefaction similarity solutions for S^, H{±), or others 



As a natural next step, before proposing concepts on existence, uniqueness, and en- 
tropy description of solutions, one needs to get a detailed understanding of the types of 
singularities that can be generated by NDEs. As often happens in nonlinear evolution 
PDEs, the refined structure of such bounded or unbounded shocks can be described by 
similarity solutions, which thus we have to begin with. 

3.1. Finite time blow-up formation of the shock wave S-{x). To this end, we use 
the following similarity solution of the NDE (11.11) : 

(3.1) u_{x,t) = g{z), z = x/{-t)'3, 
where g solves the ODE problem 

(3.2) {gg')" = lg'z in M, /(^Oo) = ±1. 

By translation, the blow-up time in (13. ip reduces to T = 0, so that we want, in the sense 
of distributions or in Ll^^ (see Proposition 13.21 below) . 

(3.3) u^{x,t) S^{x) as t — 0~. 
In view of the symmetry of the ODE (13. 2p . 



(3.4) 



9 "-^ -9, 

z —z. 



it suffices to study the odd solutions for 2; < with the anti-symmetry conditions, 

(3.5) (7(0) = /(O) = 0. 

A typical structure of this shock similarity profile gi^z) solving (13.20 is shown in Figured] 
by the bold line. 

On regularization in numerical methods. Since we are going to use essentially the nu- 
merical results now and later on, especially for higher-order NDEs, let us describe the 
peculiarities of such numerics. For the third-order equations such as (13. 2p . this and fur- 
ther numerical constructions are performed by MatLab with the standard ode45 solver 
therein. Currently, we use the enhanced relative and absolute tolerances 

(3.6) Tols = 10"^^ 

Instead of the degenerate ODE (13. 2p . we have solved the regularized equation (this idea 
is naturally associated with the general concepts of parabolic ^-approximation [TTl § 7]) 

(3.7) g'" = ^'^"^ ^ {^9'^ ~ '^9'9"^i with the regularization parameter u = 10~^^, 

where the choice of small u is coherent with the tolerances in (13. 6p . 
On asymptotics, existence, and uniqueness. One can see that g{z) is always oscillatory 
about constant equilibria =f1 as 2; ±00. Indeed, for z ^ — 1, where g{z) ^ 1, the ODE 
(13. 2p asymptotically reduces to the linear equation 

(3.8) g'" = lg'z + .... 
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m=1 : basic shock and rarefaction similarity profiles for data S^(x) 




Figure 1. The shock similarity profile as the unique solution of the problem (13. 2p . 
Hence, g{z) satisfies the asymptotics of the classic Airy function: as z — oo, 

(3.9) g{z) I Jf- cAi{z) 1 + c\z\~^ cos (^ao\z\^ + Co) , where ao = | -\/3. 

Concerning the original PDE (11. ip . this means that, as x — * —oo, and hence u —>■ 1, the 
NDE is asymptotically transformed into the linear dispersion equation 

(3.10) Ut = Uxxx, with the fundamental solution 

(3.11) 6(x,t) =r3F(a;/ti), where F = Ai(z), F" + |Fz = 0, J F = I. 

In particular, the asymptotics (13.91) implies that the total variation (TV) of such so- 
lutions of (II. ip (and hence u_{x,t) for any t < 0) is infinite. Setting = v in the 
integral yields 

+00 oo 

(3.12) |^7(-)Itv= / \g'{z)\dzr~.ji^dv = oo. 

— oo 

This is in striking contrast with the case of conservation laws (I1.17p . where finite total 
variation approaches and Helly's second theorem (compact embedding of sets of bounded 
functions of bounded total variations into L°°) are key, [26]. In view of the presented 
properties of the similarity profile g{z), the convergence in (13. 3p takes place for any a; G M, 
uniformly in R \ (5, 5), 5 > small, and in Lf^j,(R) for p e [1, oo); see below. 

Before passing to more accurate mathematical treatment of the profile g{z), we present 
its regular asymptotic expansion near the origin: for any C < 0, there exists a unique 
solution of the ODE ([S2D, (IS3D satisfying 

(3.13) g{z) = Cz + ^z^ + ^^z^ + ... (C = -0.51... for see Figure [I]). 
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m=1 : various odd shock similarity profiles g(z) for z<0 
3 I 1 1 1 1 < — 




Z 

Figure 2. The odd shock similarity profiles g{z) for z < with various limits 
g{z) C- as z — > — oo. 

The uniqueness of such asymptotics is traced out by using Banach's Contraction Principle 
apphed to the equivalent integral equation in the metric of C{—6, 6), with S > small. 

In addition, we will use the following scaling invariance of the ODE in (13.21) : if gi{z) is 
a solution, then 

(3.14) daiz) = a^(7i(^) is a solution for any a ^ 0. 

Proposition 3.1. The problem {\3.^ admits the unique shock wave profile g{z), which is 
an odd analytic function and is strictly positive for z < 0. 

Uniqueness follows from the asymptotics (13.131) and scaling invariance (I3.14p . Global 
existence as infinite extension of the unique solution from z = 0~ follows from the ODE 
(13. 2p . which, besides solutions with finite limits as z ^ —oo, admits the unbounded 
solutions with the behaviour 

(3.15) g{z) = ^z^ + ... for j^] > 1. 

Figure [2] shows other profiles from this one-parameter family with different limits of g{z) 
a.s z —oo. Of course, all of them are related to each other by scaling (I3.14p . 

The proof of the positivity is rather technical and will be focused on later. In partic- 
ular, our reliable numerics with enhanced accuracy show the exhaustive positivity of the 
similarity profile g{z), and, actually, this can be treated as a computational-based proof. 

Thus, we need to prove analyticity at z = only, since by positivity g{z) is analytic 
at any other point. Differentiating (13.21) k — 2 times yields by Leibnitz formula that the 
derivatives at z = satisfy the recursion relation 

{k + l)g'g^^) = 2{k - 2)^?('=-2) _ ^ ^W^(fc+i-.)_|±i^. 

i=2 
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It is not difficult to conclude that the growth rate of the derivatives |5'*-'^-'(0)| is not essen- 
tially more than k\. This means that the corresponding power series 

(k) 

has the unit radius of convergence. 

We now return to the convergence (13. 3p . 

Proposition 3.2. For the shock similarity profile g{z) from Proposition \3.1[ i\3.3i\ : 

(i) does not hold in L^(R), and 

(ii) does hold in Ll^^( 



Proof, (i) follows from (13.121) . since g{z) — 1 ^ L^(R_). (ii) Here we need an extra 
estimate by using the asymptotics (13. 9p : for a fixed finite / > 0, as t — >■ 0~, 



f \g{z) -l\dx = {-ty/^ J \g{z)-l\dz 

-I -«(-t)-l/3 



(3.16) ~ (-t)i/3 J |2;|^3|cos(ao|z|^)|d2: 

-/(-t)-l/3 
/3/2(-t)-l/2 

^ {-ty^'^ J v-^cos{aov)\dv {-t)T2 0. □ 



Thus, the rate of convergence in (13.161) and the fact of convergence itself depend on 
the delicate asymptotics of the Airy function, and more precisely, on the tail structure 
of the fundamental solution of the corresponding linearized equation (13.10p . Therefore, 
for NDEs such as (II. ip or other higher-order and non-fully divergent ones, the topology 
of convergence cannot be obtained in a unified manner and is individual. For instance, 
this can be L^^^-convergence, where p > 1 may depend on the NDE under consideration. 
Fortunately, due to (13.160 (note that the "gap" for such convergence expressed by the 
positive exponent ^ is rather small), for (II. 10 . we still have the Lj'^^^-convergence as for 
the conservation law (I1.17p . where this topology is naturally reinforced by Helly's second 
theorem on compact embedding of BV into 

3.2. On formation of non-symmetric final time profile. This analysis is harder but 
repeats the above arguments. Figure [3] shows a few of such similarity profiles g{z) that, 
according to (13. ip . generate non-symmetric step-like functions, so that, as t — >• 0~, 



[ C_ > for X < 0, 

(3.17) u_{x,t)^<Co forx = 0, 

( C+ < for a; > 0, 

where C_ ^ —C+ and Cq 7^ 0. In order to understand the whole variety of such non- 
symmetric profiles, one needs to check how many regular orbits can pass through any 
singular point z = zq > 0, at which g{zo) = (2:0 = 5 in Figure [3]). By Banach's 
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m=1 : shock similarity profiles with g{z^)=0, g"(Z|^)=Zjj/9, z^=5 




(gg')"=g'2/3 



-60 -50 -40 -30 -20 -10 



10 20 30 



Figure 3. Non-symmetric shock similarity profiles g{z) with various limits 
g{z) C± as z — > ±00. 

Contraction Principle, it can be shown that, for any zq > 0, there exists a ID family (a 
bundle), with the regular expansion 

(3.18) g{z) = C{z-zo) + fl{z-Zof + j^{z-zoy^ + ..., where C < 0. 

For zq = 0, this coincides with (13.131) . 

Thus, the total family of profiles regularly (moreover, analytically) passing through 
singular points z = zq is 2D, with parameters zq and C. By scaling invariance (I3.14p . 
this variety can be reduced to a ID manifold with, say, fixed zq = 1 and arbitrary first 
derivative C < 0. 



3.3. Shock similarity profiles with finite interfaces: weak discontinuities. These 
correspond to C+ = in (I3.17p . Especially, we are interested in checking how to get in 
the blow-up limit as t — > 0~ the reflected Heaviside function (ll.25p . The value at the 
origin H{0) G (0, 1) is obtained via integration. It follows from ( ]3.18p with C = that 
such regular weak solutions of the equation in (13. 2p have the following asymptotics: 

(3.19) g{z) = § (2; - zof + ... as z , 

where zq > is the only free parameter. Proof of existence of such proflles, which are 
positive for all z < zq, is similar to that of Proposition 13. II 

Using the additional smoothness of expansion (I3.19p . we set g{z) = for 2; > zq, thus 
creating an admissible type of proper weak discontinuity a.t z = Zq for the ODE (and the 
NDE), where the flux {gg'y{z) is continuous and g'{z) is Lipschitz. Several proflles from 
this family are shown in Figure HI The bold line corresponds to formation of the Heaviside 
shock (ll.25p . A mathematical justiflcation of the consistency of this cutting-off procedure 
(actually meaning flnite propagation of the right-hand interface) will be performed later. 
Here we note the following approximation property: the proflles with the flnite interface 
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m=1 : shooting shock similarity profiles with interface for data H(-x) 




Figure 4. The shock similarity profiles g{z) with finite interfaces at z = zq > 0. 



at a singular point z = zq 



(3.20) 



can be obtained as the limit as C+ 
of smooth solutions (as in Figure [3]). 



Note that the family with C+ < are composed from analytic solutions, though these 
have singular points where g = 0. On the other hand, we expect that g{z) with interface 
at 2;o > can be obtained as the limit /i — 0"*" of strictly positive solutions {g^{z)}, where 
gfj.{z) ^ /i > as z — *■ +oo. This demands extra asymptotic analysis of matching the 
flow (13.21) about g = +1 as z —oo with the ID stable manifold of (7 = as z ^ +00. 

As an alternative. Figure [5] shows another way of approximation of g{z) with zero at 
z = Zohy a. family of strictly positive analytic solutions of the ODE (13.21) . This convinces 
that the solutions with the singular point a.t z = zq and behaviour (13.191) (but not with 
finite interface!) are G-admissible in the ODE sense; see Section O 



-u, 
-t. 



3.4. Rarefaction similarity solutions. Using the reflection symmetry 
(3.21) 

of the PDF (11.11) . we have that it admits similarity solutions defined for all t > 0, 

(3.22) u+{x,t) = g{z), with z = x/t^, 
where now g{z) solves 

(3.23) {ggj' = -lg'z in R, /(^oo) = Tl. 

Obviously, these profiles are obtained from the blow-up ones in (13. 2p by reflection. 



(3.24) if shock profile g{z) solves (13.21) . then g{~z) is rarefaction one in (13.231) . 
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m=1 : construction of shock similarity profile by analytic approximation 




Figure 5. Analytic approximation of g{z) with the behaviour ()3.19p at z = 
ZQ > 0. 

The corresponding rarefaction similarity profile is shown in Figure [H This regular simi- 
larity solution of (13.231) has the necessary initial data: by Proposition I3.2( ii). in Lj^^, 

(3.25) u+{x,t) ^ S+{x) as t-^0+. 

Other profiles g{—z) from shock wave similarity patterns generate further rarefaction 
solutions including those with finite left-hand interfaces. 

Thus, all three Riemann's problems for the NDE (II. ip with data (ll.20p and (11.250 
admit a unified similarity treatment. 

4. Unbounded shocks and other self-similar singularities 

We continue to introduce other types of shocks and singularities that are associated 
with the NDE ffTTIl . 

4.1. Blow-up self-similar solutions: invariant subspace and critical exponent 
= —jq- We now consider more general blow-up similarity solutions of (II. ip . 

(4.1) u^{x,t) = i-trgiz), z = x/{-tY, P = ^ (aGR), 
where g solves the ODE 

(4.2) A{g) = {gg')" =^-^g'z~ag = Cg in R. 

In (14. ip . we introduce an extra arbitrary parameter a G M. We next show that the 
behaviour of the similarity profiles g{z) and hence of the corresponding solutions Ua{x,t) 
essentially depend on whether a > or a < 0. 

We first prove the following auxiliary result explaining a key feature of the ODE (14. 2p . 
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Proposition 4.1. (i) Both nonlinear A and linear C operators in [4-2) admit the 4D 
linear space 

(4.3) W^ = Span{l,z,z'^,z^}. 

(ii) The ODE 1 14 -S^ possesses nontrivial solutions on W4 in two cases: 



(4.4) (I) a = ac = —jq, with the solutions given by 



(4.5) g{z) = Cq + CiZ + ^ z , where Co,i G M are arbitrary constants, and 

(4.6) (II) a = -l, with g{z) = ^Cl + 20Clz + C2Z^ + ^z\ C2eR. 
Proof, (i) is straightforward, since, for any 

(4.7) g = Co + Ciz + C2Z^ + C^z^ e W^, 

the following holds: 

A(^) = 6(CiC2 + C0C3) + 12(C| + 2CiC3)z + GOCaCg^^ + 60Clz^ G W4, 
^^■^^ Cg = -aCo + ^Ciz + ^^ Caz^ + Csx^ e W4. 



(ii) According to the ODE (14. 2p . equating the coefficients given in (14.81) yields the 
algebraic system 



(4.9) 



6(^1^2 + C0C3) = -aCo, 
60C| = C3. 



The last equation yields C's = ^ (we exclude the easy case C3 = 0), and then the third 
one implies that either C2 = 0, or 

(4.10) 60C2^ = ^C2 =^ a = -l. 

Assuming first that C2 = and substituting into the second equation in (14.91) . we infer 
that, for arbitrary Ci 7^ 0, 

(4.11) 2AC,i, = ^ I = i.e., a = 



The first equation is also valid for any Cq € M. 

Choosing C2 7^ leads to the less interesting case (14. 6p . □ 

Remark: blow-up for the NDE on the invariant subspace. Since the linear 
subspace (14.31) . W4 = Span{l, x, x^, x'^} is invariant under the quadratic operator A, 
the whole original NDE (11.11) can be restricted to W4, on which the evolution of the 
corresponding solutions, 

(4.12) u{x, t) = Co{t) + Ci{t)x + C2{t)x^ + C3{t)x^ G W4, 
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m=1 : shock similarity profiles for various a>0 and a<0 



a=0.1 



g'(0)=-o.5i 




(gg')"=(1+a)g'z/3-ag 
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Figure 6. The shock similarity profiles as solutions of the ODE (14. 2p for a > 
and a < 0. 

is governed by a dynamical system with the right-hand sides from (I4.9p . 



(4.13) 



C'o 



6(CiC2 + CoC3), 

12(C| + 2CiC3), 
6OC2C3, 



This is easily integrated starting from bottom and gives the blow-up behaviour, 



60{T-t) ' 



Do 



+ 



(T-t)2/5 
72OA0 



(T-t)l/10 ^ (T-t)2/5 (T-t)3/5' 

where the blow-up time T, Aq, -Bq, and Co are arbitrary constants. Unlike the above 
gradient blow-up, the present singularity formation is governed by the cubic growth of 
the solution fl4.12p as x — > 00, so that it is less interesting. 

4.2. Blow-up similarity solutions and their properties for a G (ac,0). Basic local 
mathematical properties of the ODE (14. 2 p are similar to those for (13. 2p . so we omit these 
details and concentrate on principal global features of blow-up patterns. The typical 
structure of solutions g{z) of (14. 2p can be understood from Figure O where we present a 
few profiles for a > and a < 0, including the bold line for a = that has been studied 
in the previous section. All these profiles satisfy the anti-symmetry conditions (13.50 . 

For a > 0, the ODE (14. 2 p suggests (and Figure [6] confirms this) that the similarity 
profiles are growing as z — > —00 according to the linear part of the equation, i.e.. 



(4.14) 



for 2: < -1 



3a 

\z\ 1+" 



+ 



C„ > 0. 
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Passing to the limit t ^ in (14. 11) does not lead to a shock wave but to a "weak 
singularity" , 

3a 

(4.15) Ua{x,0^) = C^\x\~ for X < and a; > 0. 

For a = I or a = 2, these give the behaviour ~ x or ~ near the origin respectively. 

The behaviour for a < is more interesting, and deserves further analysis. First of 
all, for a G (— ^^,0), the similarity profiles g{z) are still positive for 2; < and have 
the asymptotic behaviour f l4.14p . Therefore, (14.151) is also true but one should take into 
account that here the exponent < 0, so that the final time profile is unbounded, i.e., 

(4.16) Ua{x, 0^) = C^\x\^ ^ ±00 asx^O'F. 

This is the first example of strongly discontinuous unbounded shocks that can be obtained 
evolutionary via blow-up of sufficiently smooth continuous self-similar solutions of the 
NDE (11.11) . The conservation law (11.171) cannot produce such a shock, since by the 
Maximum Principle, 

sup^, \u{x,t)\ < sup 2; |'Uo(2;)| for all t > 0. 

In Figure [71 we show how the positive similarity profile g{z) for 2; < is deformed 
according to the limit 

(4.17) g{z) g^{z) as a ^ a' = 

Note that a slightly easier mechanism of formation of such "saw cusps" is a well-known 
bifurcation phenomenon in the FFCH (11.331) and Degasperis-Procesi (11.361) type equa- 
tions; see [371 P- 422] and references therein. Figure (Ht^a) explains the actual convergence 
of positive profiles in the limit (I4.17p . while (b) shows the oscillatory structure of gc{z) 
for z <^ —1 a = 

Finally, in Figure [9l we show other examples of non-symmetric profiles g{z) for various 
a including those in (b) that have the finite interface fixed at Zq = 2.192. The bold profile 
generates as t ^ 0~ the Heaviside function (11.251) . with parameters 

(4.18) //(O) = 0.4197..., 20 = 2.192.... 

Remark: complete blow-up. There is an important question concerning continua- 
tion of the similarity solution satisfying (I4.16P beyond blow-up, i.e., for t > 0. This is 
not the subject of the study, so we present a few comments on that. According to ex- 
tended semigroup theory of blow-up solutions (see [9] for references and known results 
for reaction-diffusion PDEs), it is natural to construct such an extension via a smooth 
approximation (truncation) of the NDE by introducing a bounded nonlinearity, 

(4.19) M„ : Ut = {^n{u))xxx, with v'„(u) = i | as n ^ 00. 

Taking the same initial data we obtain a sequence {un{x, t)} of better bounded solutions 
(bounded since blow-up demands the quadratic growth ~ as m — > 00 of the coefficient 
that is not available anymore; we omit difficult technical details). 
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m=1 : shock similarity profiles fora=0, -0.02, -0.05, and -0.09999 




the "saw" profile). 



m=1 : Formation of oscillatory profile as a tends to a =-1 /1 m=1 : local structure of oscillations for z — 1 00, a =-1 /1 

, , , , , , 1 0.4 1 I , , , 




(a) convergence to gc{z) (b) zeros of gc{z) 

Figure 8. Tlie limit (|4.17p : convergence to tlie singular profile gdz) (a), and 
the asymptotic structure of the "tail" of gdz) for z ~ —100. 



Nevertheless, each Un{x,t) can produce a shock at some t = t„ ~ 0, but now of finite 
height ~ 2n. Since Unt{x,tn) gets very large in a neighbourhood of x = 0~ (this follows 
from the fact that, for the original solution, Ut{x, 0^) +oo as x 0^ for "data" fl4.16p ). 
one can expect that Un{x,t) for x < grows rapidly for small t > achieving fast the 
value ~ n. This means that Un{x,t) for small t > gets the shape of nS^{x). 
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m=1 : formation of singular profile as a tends to -1/10 from above 



m=1 : similarity profiles with finite interfaceat z =2.192 for various a 
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(a) C_ 7^ —C+ (b) with finite interface 

Figure 9. Various similarity solutions of (j4.2p for positive and negative a. 

Finally, we then conclude that for any fixed t > 0, as n — > oo, 

-oo for X < 0, 
■oo for X > 0. 



(4.20) 



t) ~ nS^{x) 



Actually, this means that this function u{x,t) = ±oo is the correct proper solution as a 
natural continuation of the similarity solution (14. ip . (14.160 for t > 0. In other words, we 
observe complete blow-up at t = 0~ and the fact that the extended semigroup of proper 
(minimal) solutions is discontinuous at this blow-up time. 

4.3. Oscillatory invariant profile for ac = —Jq- Then the ODE (14. 2 p takes the form 

(4.21) {99r = lo9'^ + ro9. 

which is associated with the invariant subspace W4 in (14.31) . 

We have already seen, that the critical case a = ac is of special interest. Figure [TOlfa) 
shows that, in fact, for = — ^, there exist two continuous solutions. The second one 
given by the dotted line is the exact invariant one 

(4.22) g*{z) = -mz + ^ z'^ e W4, where - m = g'{0) < 0. 

Figure [TUT b) explains the structure of the envelope to the main "saw-type" oscillatory 
profile gc{z) that has the same structure as in (14.140 . i.e., 

(4.23) L{z) = C-\z\-^ as z ^ -00. 

Therefore, the final time profile Ua{x,0~) is not like (14.161) . Indeed, the envelope L{x) 
only has this form, so that we have another topology of convergence: 

(4.24) Ua{x,t) ^ C-L{x) as t ^ 0" 

in the weak sense in L°° or in L^. Here C_ = 1.67... is a constant that characterizes the 
oscillatory part of this shock wave distribution; see more details below. 
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m=1 : two shock similarity profiles for a =-1/10 



m=1: Envelope to "incomplete blow-up" shock profile for a =-1/10 





(a) two profiles (b) envelope to the saw 

Figure 10. Two solutions of (fOT]) on W4. 

m=1 : a tends to a|,=-1 /1 0, formation of "incomplete" singularity -|z-Zjj| 




Figure 11. Local convergence as a — > — ^ of smooth solutions g{z) to the 
"saw" gc{z) with the behaviour (j4.25p . 

Figure [TT] explains formation of zero at some z = zo < of gc{z) as a ^ 
clearly seen that g^z) has a symmetric behaviour for z zq, so instead of fl3.18p . 

(4.25) g,{z) = C\z -z^\ + %{z- z,f + ... , C > 0. 

Since the "flux" {gl)" is continuous at 2; = zq, fl4.25p represents true weak solutions of 
fl4.2ip . Moreover, it is also crucial for our further analysis that 

the weak piece- wise smooth solution gc{z) can be obtained 
as the limit of solutions of the family of ODEs (14.20 . 



(4.26) 



Finally, we prove the following result concerning the invariant properties of the oscilla- 
tory profile gc{z): 
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Proposition 4.2. For a = —jq, in the piece-wise sense, 

(4.27) g,{z) e W^. 

Proof. Without loss of generality, we consider the solution gdz) of (14.211) with m = 1, 

(4.28) g{0)=g"{0) = and ^?'(0) = -1. 

By uniqueness, it follows from (14. 5 p that the first hump of the profile shown in Figure 
M^a) is explicitly given by 

(4.29) g^{z) = -z + j^z^ for zE {zo, 0), where Zq = -VOO {gdzo) = 0). 

For the second hump with the transition given by fl4.25p . we take the "full" profile (14. 5p 
and demand two matching conditions at 2; = 2:0, 

g{zo) = 0, g'{z^) = 2 = -giz^), i.e.. 



Co + CiZo + ^ ^0 — 0) 
^1 + M ^0 = ~2. 



(4.30) 

20 

This yields Ci = —5 and Cq = — 4-\/60, so that the second hump is given by 

(4.31) g,{z) = -AVm - 5z + ^ z^ on {zi,zo). 
The second zero zi is 

(4.32) zi = -pVeO =^ p^-5p + 4: = 0, i.e., p = = 1.56155... . 

This procedure is continued without bound and gives the "saw" as in Figure [TOT b). □ 

4.4. Local complete blow-up in the ODE for a < ac- It turns out that for a < —jq, 
the orbits g{z) cannot pass through singularity at {g = 0}. This local "complete blow-up" 
is clearly illustrated by the equation for the second derivative 

(4.33) v = g" =^ gv" + + Ag'v' = ^vz + ^g'. 

In other words, this blow-up v —00 as z ^ Zq < is performed according to the 
quadratic source term, 

(4.34) v" = -lv^ + .... 

The ODE (14. 2 p for a ^ —1 admits the following behaviour at such complete blow-up 
zeros: 

(4.35) g{z)=C./F^+^^zo{z-zof + ... as z ^ z^ , C > 0. 

For a = —1, the remainder is ~ 0{{z — 2;o)^)- The results of such shootings are presented 
in Figure [121 
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m=1 : a tends to a from below 




10- 



Figure 12. Complete blow-up of orbits for a < 
4.5. Rarefaction similarity solutions. These are defined as follows: 

(4.36) Ua{x,t) = t"g{z), with z = x/t~, where 

(4.37) {gg'y = -'-^9'z + ag in R. 

This gives the reflected shock blow-up profiles g{—z), and all the above properties can be 
translated to the rarefaction evolution. 

In connection with this, we mention that the unbounded initial data (see (14.161) ) 



(4.38) 




for a; > 0, 
for X < 0, 



for any a G (— ^Ij^O) generate a bounded (for any t > 0) similarity solution (I4.36P with 
the corresponding profile g{—z). According to our further theory, this is an example of 
an entropy solution that is expected to be unique. 

5. On G-admissible solutions of ODEs 



The approximating properties that we underlined in (I3.20p and (14.261) are key for under- 
standing the ODEs involved. Actually, these are related to the concept of G-admissibility 
of solutions of ODEs with shocks. It was introduced in 1959 by Gel'fand [TDj and was 
developed on the basis of TW-solutions of hyperbolic equations and systems; see details 
in [m § 2, 8]. 

Namely, talking for definiteness the typical ODEs (13. 2p or (14.20 . a non-classical (i.e., not 
C^) solution is called G-admissible if it can be constructed by a converging sequence of 
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smooth solutions of the same ODE. In particular, this means that the profiles with finite 
interfaces as in Figure [Hand even the "saw-type" profile Qc for G-admissible. 

Of course, this is an ODE concept, but we will bear it in mind when develop a PDE 
concept of entropy solutions in [131 IE] • 



6. On travelling wave and generic formation of moving shocks 

We now briefly discuss how to construct moving shock waves, since previously we 
concentrated on standing shocks such as S'_(x), H{—x), and others. Of course, the first 
simple idea is to consider: 

6.1. Travelling waves (TWs). These are solutions of (11.11) of the form 

(6.1) u{x, t) = f{y), y = X- Xt, 

where A G M is the TW speed. Substituting into (II. ip and integrating once yields 

(6.2) -Xf = {ff)" =^ _xf = {ffJ + Ao, 

where G M is an arbitrary constant. If a shock (a finite discontinuity) occurs at some 
point yo G M, (16.21) yields the following Rankine-Hugoniot condition for the speed: 

(6-3) A = -Mmfl, 

where, as usual, [(Ollz/o) stands for the jump of the function (■)(?/) at y^. In the case 

(6.4) A = = 0, 
we obtain either constant solutions 

(6.5) fiy) = C, and hence / = S±{y), H{±y), and others, or 



(6.6) f{y) = ,/Aiy + A2, where Ai,2 G M, 

from which one can reconstruct various discontinuous solutions satisfying the Rankine- 
Hugoniot condition (16.31) . At yo = (cf. (I4.35P ). the solution (16.60 is not sufficiently 
smooth to be treated as a solutions of the Cauchy problem. 
For A 7^ 0, (16. 2p admits smoother explicit solutions 

(6.7) f(y) = -l(y + By-l^Ao, Bern, 

which also give various shocks. In general, it is easy to describe all the solutions admitted 
by the standard second-order equation (16. 2p that we are not doing here and will focus on 
the following principle question: 
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6.2. On formal theory of generic formation of finite moving shocks. Let a uni- 
formly bounded continuous weak solution u{x,t) of the NDE (11.11) in R x (—1,0) (for 
convenience, the focusing blow-up time is again T = 0) create a discontinuous shock at 

(6.8) X = 0, t = 0~~ , at which the shock moves with the speed A. 

Then it is natural to consider this phenomenon within the moving frame and introduce 
the corresponding TW variable y = x — Xt as in (16. ip . Hence, the solution u = u{y,t) 
now solves 

(6.9) Ut = {UUy)yy + XUy ^/GM, ^^( — 1,0). 

We now study formation of bounded shocks associated with the similarity variable z 
in (13. ip . Indeed, there are other blow-up scenarios, with variables as in (14. ip : see [14J for 
more generic single point "gradient catastrophes" (these require other rescalings). Hence, 
according to (13. ip . we introduce the rescaled variable 

(6.10) u{y,t) = v{z,t), where z = x/{—t)3^ r = — ln(— t) ^ +oo as t ^ 0^. 
Substituting this into (16.91) yields that v{z, r) solves the following non-autonomous PDE: 

(6.11) Vt = {vvz)zz — ^VzZ + Xe'^'^gz in R x ]R_|_. 

Here the right-hand side is an exponentially small for t —>■ oo perturbation of the station- 
ary operator of the ODE in (13. 2p . Therefore, it is natural to expect that the asymptotic 
behaviour as r — -\-oo (i.e., t ^ 0^) is described by the non-perturbed equation, 

(6.12) Vr = {vvz)zz -\vzz = B(t;) = K{v) - Cv in R x M+. 
Therefore, stabilization in (16. lip and (I6.12p to a self-similar profile, 

(6.13) v{z,t) —>■ g{z) as r ^ oo uniformly on compact subsets, 

is equivalent to the question of the stable manifold of the stationary solution g{z). In the 
linear approximation, this will depend on the spectrum of the non self-adjoint third-order 
operator 

(6.14) B'ig)Y = {gYy"-lY'z. 

This spectral problem is not easy, though similar higher-order linear operators already 
occurred in some evolution odd-order PDEs; see [Hj § 5.2], [ini § 9.2] for details. 

Considering, for simplicity, B'{g) in the space of odd functions defined in R_ = {z < 0}, 
in view of conditions (13. 5p on g{z) and the smooth behaviour (I3.13p . the end point z = 
can be treated as a regular one. At infinity, where g{—oo) = 1, we arrive asymptotically 
at the linear operator as in (13. 8p . which admits Airy function as one of its eigenfunctions. 
These operators are naturally defined in the weighted space 

(6.15) Ll{R+), where p(;z) = e""!"!''" 

and a > is a sufficiently small constant; see more details in flU^^ § 9.2]. 
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It is important that, in view of the scahng symmetry (13.141) . 

A = G cr(B'(gf)) with the eigenf unction 

(6.16) 

so that B'(5f) has the non-empty centre subspace E'^ = Span{?/^o}- 

Then, according to typical trends of invariant manifold theory; see [25] (we must admit 
that sectorial properties of B'{g) in are still unknown, while compactness of the resol- 
vent is most plausible), it is natural first to try a centre subspace behaviour for (I6.12p . 
Stability on the centre manifold is key for understanding the stabilization properties of 
the flow (I6.12p . To this end, we perform the "linearization" procedure 

(6.17) v{z,T)=g{z)+Y{z,T) 

to get the following equation with a quadratic perturbation: 

(6.18) Yr = B'{g)Y + l{Y^),,,. 

In order to formally check the stability of the centre manifold, we assume that for r ^ 1 
the behaviour follows the centre subspace, i.e., 

(6.19) Y{z,t) = a{T)ipo{z) + w{z,t), where w.LE'^. 

Here, the orthogonality is defined according to the corresponding Eucledian space with 
the indefinite metric (■,■)*, 

{v,w)^: = Jv{z)w{—z) dz; 

see details and examples in [101 § 9.2]. Substituting (16.191) into (I6.18P and multiplying by 
the corresponding normalized adjoint eigenf unction ip^ yields 

(6.20) a' = -700^ + ... as r ^ cx), where 7o = |((^o)"'> V'o)*; 

where the scalar products involved are defined by the extension of the linear functionals 
by Hahn-Banach's classic theorem (we omit the details that are not crucial here). For 
existence of such a centre subspace behaviour, it is necessary that 

(6.21) 70 ^ =^ 3 a(r) = -j- + ...^0 as r ^ oo. 

On the other hand, stabilization in (16.131) is also possible along a stable infinite- 
dimensional subspace of B'{g)- Note that this would be the only case if 70 = in 
(16.211) . Such operators (without degeneracy at z = 0, which nevertheless is not strong 
and/or essential) are known to admit a point spectrum [TOl § 9.2]. Such formal stable 
manifold behaviours admit similar standard calculus, though the duality between the op- 
erator B'{g) and the adjoint one B'*{g), as usual for odd-order singular operators, will 
provide us with several technical difficulties. Recall again that there exist other more 
structurally stable (and possibly generic) ways to create single point gradient blow-up 
singularities, where other rescaled nonlinear and linearized operators occur, [T^. We do 
not touch those questions here. 

Overall, this explains some formalities of the asymptotic stability analysis, which re- 
mains a difficult open problem. Even checking the necessary non-orthogonality assump- 
tions, such as (16.211) . can be a difficult problem, both analytically and numerically. 
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Anyway, it is key that if stability of stationary profiles g{z) for (16.121) and flG.lip takes 
place (at least, on some manifold) then, for bounded moving shocks, we can use a similar 
geometric rule (11.221) distinguishing 5-entropy and non-entropy discontinuous solutions 
(see more details in [12]). 
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